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We consider the physics of a matrix model describing DO-brane dynamics in the presence of an 
RR flux background. Non-commuting spaces arise as generic soltions to this matrix model, among 
which fuzzy spheres have been studied extensively as static solutions at finite N. The existence of 
topologicaly distinct static configurations suggests the possibility of D-brane topology change within 
this model, however a dynamical solution interpolating between topologies is still somewhat elusive. 
In this paper, we study this model in the limit of infinite dimensional matrices, where new solutions- 
the fuzzy cylinder and the fuzzy plane among them- appear. We argue that any dynamics which 
involves topology change will likely only occur in this limit, after which we study the decay of a fuzzy 
cylinder into an infinite collection of fuzzy spheres as both a classical and a quantum phenomenon. 
We conclude from this excercise that in certain limits, matrix models offer a viable framework in 
which to study topological dynamics of fuzzy spaces, and could perhaps be a precursor to a viable 
theory of space-time topological dynamics. 

<: 

0> 1 I. INTRODUCTION AND MOTIVATION 

[ Non commutative geometries have long been studied extensively in their own right as a possible framework 
■ for spacetime physics at very high energy scales 29] . From their genesis as an ad- hoc manner in which to at- 
tempt a regularization of quantum field theory non commutative geometries have since been motivated by 
semi-classical studies in general relativity 8] and string theory. As an area of intrinsic mathematical interest [l| 
it apparently needs no justification. More recently such geometries have found a natural realization as solutions 
of various matrix models describing D-brane dynamics, (see 3@E3 an< ^ refernces therein). It is the aim of 
this report to study a particular matrix model describing DO-brane configurations in the presence of RR flux 
as a concrete example in which one finds interesting time dependent physics for these non commutative spaces, 
including possible topology changing transitions between distinct geometries. By topology in the context of non 
commutative spaces, we specifically refer to the properties of the commutative manifold associated with the fuzzy 
geometry we are considering. Hence by topology changing dynamics, we are refering to dynamics which in the 
commutative limit (if such a limit is available), would look like topology change in the classical sense. It appears 
that this putative aspect of matrix dynamics is closely related to its nature as a theory with U(N) invariance. 
The eventual hope is that matrix models could provide a realistic framework in which to study the toplogical 
dynamics of extended objects embedded in spacetime, if not spacetime itself- a widely anticipated (if sometimes 
controversial) corollary of any putative quantum theory of gravity, and one that is deserving of a brief review. 
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Ever since Wheeler's initial musings on the matter 0, the notion of fluctuating spacetime topology lay more 
or less dormant in the minds of theoretical physicists for several decades. Although most physicists have long 
agreed that notions of spacetime would have to be radically reformulated at the Planck scale, we have only just 
begun to uncover the physics likely to be responsible for this reformulation. Inspite the multitude of approaches 
to modelling spacetime at the Planck scale (among them, canonical and covariant quantum gravity, loop quantum 
gravity, string theory and non-commutative geometry), common to them all is the notion of a minimally resolvable 
length scale. This notion is either taken up as an a priori input to the theory, or is a derived consequence of other 
principles. For instance, in loop quantum gravity, it is a non-trivial result of this theory that various geometric 
operators such as volume and area are quantized in units of h 19] 30] . In string theory, T-duality, the superstring 
uncertainty principle and the appearance of non-commutative geometry (NCG) in certain situations all imply a 
minimal length scale, though its interpretation depends on precisely which of the three mechanisms are at work 31] . 
In NCG, one takes as a starting point the non-commuting nature of space (and/or time), which is logically equivalent 
to positing a spacetime uncertainty principle and formulates physics afresh in this framework. Non-commutative 
geometry naturally makes an appearance in string theory through the dynamics of D-branes |3]|5l]|l8j] | l3 | and 
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through the theory of relativistic membranes It will almost certainly have a central role in the future if 

the 'M' in M-theory turns out to stand for matrix. Hence from the perspective of a string theorist, NCG would 
be the most natural and best motivated avenue in which to directly study the nature of spacetime at the Planck length. 

After the hard work of consistently formulating any putative model of spacetime at the Planck length has been done, 
the question of whether or not topology changing dynamics is a feature of this model, or whether the theory breaks 
down before such dynamics can take place, immediately beckons. In considering this question, perhaps it is appropri- 
ate to refresh ourselves with what the presently accepted paradigm of space time dynamics (classical general relativity) 
has to say. It should be clear that general relativity (derived from the Einstein-Hilbert action), does not contain any 
scope for topological dynamics. One can do general relativity in a background of any given topology, but the resulting 
dynamics does nothing to change this topology. This is easily seen from the fact that constructively, any topology 
changing process will involve the formation of singularities, which is where general relativity breaks down. Further- 
more, it has been demonstrated [23| that general relativity excercises a sort of 'topological censorship' principle, in 
that at intermediate energies (below the compactification scale and above the Hubble scale) all causal curves (probes) 
that traverse non-trivial toplogical features of spacetime causally disconnect from the causal curves (observers) which 
do not. That is, the most we will ever observe of a non-trivial topological feature of spacetime is a singularity. 

Adopting a 'top down' approach, that is one which bases itself on established physical paradigms, euclidean quantum 
gravity posits a sum over histories approach to quantizing general relativity. Topology changing processes have been 
extensively studied in this context (see [22] for a review and references to the relevant literature). By this we mean 
that within the sum over histories, histories which interpolate between 3-d manifolds which are non-diffeomorphic 
to each other (cobordisms) are analysed for their likely contribution to the path integral. This program has had 
some degree of success, particularly in ruling out certain classes of acausal topological transitions (cobordisms) [2^ 
within this framework. However the approach still has a lot of fundamental issues to answer for, not least of which 
being the definition and consistency of the sum of histories (path integral) approach to quantizing general relativity. 

String theory on the other hand, has already conclusively demonstrated the capability to model spacetime 
topology change (^3|@| and references therein). Starting from the observation that the Kahler moduli space 
(which defines the Calabi-Yau 3-fold on which the theory is compactified, and hence the CFT describing it 
on the worldsheet) separates into distinct cells which define distinct compactifications, Aspinwall, Greene and 
Morrison [2^ showed that there was no physical obstruction to marginally deforming the underlying worldsheet 
CFT such that it moves into a different cell in moduli space. This has the effect of changing the target space of 
the theory in such a way that in certain cases, the topology of the target space has changed. Such geometrical 
transitions are continuous and are known as 'flop' transitions (previously known to algebraic geometers), and 
have been extensively studied since. However it is in keeping with the perspective natural to string theory- 
that geometry is notion derived from the dynamics of our degrees of freedom- that the study of such transi- 
tions involves the geometry only indirectly. The calculations involved are quite formidable, and although this 
remarkable result is not to be detracted from, if our primary concern is to directly model quantum spacetime 
dynamics, one would be tempted to look for another avenue to explore this question. This temptation only becomes 
stronger when one would like model more than just the dynamics of the internal space of a given compactification. 

Non-commutative geometry seems to offer exactly such an avenue, and as delineated by Madore and SaegerpH 
should rather naturally describe topology changing dynamics, once the dynamics have been determined. Matrix 
models (in a certain reference frame) naturally offer just such a dynamical framework. We stress that the particular 
framework that we are about to work in, will specifically concern examining D-brane topology change. We will be 
studying the dynamics of non commuting spaces which appear as solutions to a rather well known matrix model 
describing DO-brane dynamics, with the aim of demonstrating the existence of dynamical topology change in this 
context. One only has to realise the possibility that our universe might be a 3-brane embedded in a higher dimensional 
spacetime to appreciate the relevance of this study towards the eventual goal of understanding spacetime topology 
change. 

Henceforth, any refernce to topological dynamics will specifically refer to dynamics in the context of non- 
commutative geometry. However even within this domain, there are different approaches in the literature concerning 
the type of topological transition studied. We wish to point out that to our knowledge, the type of non-commutative 
topological transition we wish to report on is dictinct from the types hitherto looked at (some of which have demon- 
strated promising preliminary results- see for example 30 ^3 [23] ) • To understand this, it helps to first categorize 
the types of possible topological transitions thought possible in the context of non-commutative geometry. 



3 



A. Classes of Topological Transitions 

We start with the basic observation that matrix models often have as solutions to their (static) equations of motion, 
matrices that define an algebra. In the event that this algebra is a Lie algebra, they satisfy the relation: 



[Xi, Xj] = ifijkXk (1) 

That is, any set of matrices that satisfies the above is a solution to the matrix model in question. The structure 
constants fijk define the Lie algebra. In general, for any given dimension N of the matrices in question, there are 
many solutions to Q which correspond to the number of ways N can be partitioned into the integers defining the 
irreducible representations of the Lie algebra. That is, one could have any direct sum of reducible representations 
consistent with the total dimension of the matrices being equal to N. Of course, one also has the irreducible solution 
(if such an N dimension irrep exists). Furthermore, if (JJJ also describes the structure relations of a non-commutative 
geometry, one can interpret this as saying that the matrix model has solutions corresponding to multiple copies of 
the fuzzy spaces defined. To illustrate this, let us hypothetically take the index i to run from 1 to 3, and the fijk to 
be propotional to the 3-d Levi-Civita tensor: 



[Xi,Xj] = inEij k X k (2) 

In which one recognizes the X, to define the structure relations of the fuzzy sphere (see |T^ for a review of its 
construction), with the non-commutativity parametrized by k. Any representation of su(2) (irreducible or otherwise) 
given by the generators { Ji} furnishes the solution: 

Xi = nJi (3) 

For an irrep, the casimir operator XiXi — K 2 j(j + 1) gives us the square of the radius of the fuzzy two sphere, 
where the dimension of the representation is given by N = 2j + 1. Any reducible representation can be expressed as 
a direct sum of irreps: 

X % = ® rK jf ) (4) 

Where the summation r runs over all the irreps contained in the representation. This solution de- 
scribes a set of multiply superimposed spheres. The radius of each constituent 2-sphere is given by 
= n2 jr(jr + 1) and the dimension of the representation is now given by N = J2 r (^jr + 1). 

Although the equations of motion do not fix the representation, the energy of the solutions corresponding to different 
representations will differ in general. The example we are presently studying can be derived from the matrix model 
defined by the action: 

S = T J dtTrf^XiXi + ^[X^X^X^Xj] - (5) 
We will introduce this model more formally in the next section. The energy of the static solutions is given by: 

V = Tr(-^[Xi,Xj][Xi,Xj] + -KtijkXilX^XM]) (6) 
From which we see that the energies of our fuzzy sphere solutions are given by: 



E = -T K 4 -Y,Mjr + 1)(2> + 1) (7) 

r 

Where again the summation is over the irreps contained in the solution. We can immediately conclude 

that for a given N, the irreducible representations correspond to the lowest energy (bound state) solutions. 
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The observation that reducible representations have a higher energy has motivated the study of transitions 
between these representations. Geometrically, the simplest example of this would correspond to two spheres 
melding into one (^lJB an d references therein). This would certainly involve a change in topology, as two 
disconnected spaces fuse into one connected space. We shall refer to this type of topological transition as l iso- 
geometric\ That is the algebraic and geometric structure (0) of the solutions does not change in the transition. 

A second class of topological transitions was studied in 9] , where the observation that a given Lie algebra structure 
could give rise to several distinct fuzzy geometries was used to study possible transitions between these geometries. 
The example considered was the case when Q defined the Lie algebra of SU(4), which contains fuzzy versions 
of CP(3), CP (2) and S 2 x S 2 , among others Transitions between these fuzzy geometries were explored and 
obstructions to topology change were uncovered arising from the equations of motion which fixed the particular 
representation of SU(4). Nevertheless, this work does highlight the possibility of a second type of topological 
transition, which we shall classify as l isoalgebraic\ That is, that although the geometry and topology of the fuzzy 
spaces might undergo a transformation (e.g. CP (2) — ► S 2 x S 2 ), the algebraic structure of the matrix solutions JU is 
preserved in this transition. Although an isogeometric transition is also isoalgebraic, the converse need not be true. 

In the present work, we wish to study a third class of topological transitions which are neither isogeometric nor 
isoalgebraic. By this, we mean a transition such that the algebraic structure of the solutions changes in time, with 
the initial algebra describing one fuzzy geometry and the final algebra, describing another. The specific example 
we will be studying will be the decay of a fuzzy cylinder into an infinite collection of fuzzy spheres. This study 
is motivated by the observation that a rather well known matrix model which we are about to introduce, contains 
solutions corresponding to different algebras in the limit of infinite dimensional matrices [l^. This limit is neccesitated 
by the fact that some of the solutions of this model only have infinite dimensional representations. Although going 
to such a limit might seemingly involve a loss of calculational ease, for the questions we are interested in, we will find 
that this is not the case. In fact we feel that there is a case for why one is likely to encounter toplogical dynamics 
only in this limit, and in a sense this limit could be a prerequisite for our investigation of topology change. 

B. Topology change only in the Large N Limit? 

A generic feature of matrix models of interest, is their global U{N) invariance. The action defining our example 
|J5J is invariant under the symmetry Xj — > UXiW, and it is this feature alone of matrix actions that has generated 
much theoretical interest. It has been observed that as N tends to infinity, the symmetry group U(N) tends to the 
area preserving diffcomorphism groups of certain 2-dimensional surfaces, such as S 2 or T 2 . This association arises 
in the matrix regularization of membrane theory, where the original diffcomorphism invariance of the worldvolumc 
manifests as a U(N) invariance of the associated matrix model. However, it turns out that one ends up with the 
same U(N) invariance when one regularizes membranes of any given genus This is due to the fact that the group 
of area preserving diffeomorphisms of any 2-dimensional surface (of arbitrary genus) can be approximated by U(N) 
in the large N limit j^]. That is, the matrix regularized theory appears to have more structure than the theory it 
was meant to approximate by this virtue of U (N) , and seems to describe a theory where the membrane topology can 
change. From this observation we feel that it is highly plausible that any matrix model with U (N) invariance could, 
and should only describe the topological dynamics of fuzzy spaces in this large N limit. This observation seems to 
have borne out in practice as the example we are about to study, (which involves infinite dimensional matrices) has 
proven to be remarkably tractable when compared to studying similar questions in a finite dimensional setting. 

II. THE MODEL 

The matrix model alluded to previously, which we will be study ing is described by the following action, describing 
DO brane physics in the presence of an RR flux background (30[3E3E3I13l) : 

S = T J dtTr(^XiXi + ^[Xi,Xj][Xi,Xj] - l -ne ljk X l [X. J , A fe ]) (8) 

We follow the conventions of and work in units where 27ra' = 1. The Xi are taken to be N x N, hermitian 
matrices which can be taken to be traceless, and the DO brane tension is given by Tq = \/2n/g s . The index i runs 
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from 1 to 3. This action is supplemented by the Gauss law constraint: 

[X i ,X i ]=0 (9) 

Which arises from the Aq equation of motion for the RR gauge field. Note that the last term in the action is a 
Chern-Simons term which induces interactions between the DO-branes through the 4-form flux, which assumes the 
vacuum expectation value Fqvo 3 = —2ncijk- This term was deduced by Myers by demanding consistency of the D- 
brane action with T-duality |5j . This model also appears in several other contexts outside of D-brane physics (see 
and references therein). However we shall adopt the perspective of a non-commutative geometer in this report, 
and only utilise this model insofar as it provides a model for the matrix dynamics of fuzzy spaces. The equations of 
motion derived from this action are: 

X t = [X.j, ([X t ,X.j\ - iK€ ijk X k )\ (10) 

Inspite of its appearance, the above results whether or not one takes the Gauss law constraint into account (the 
only modification would be terms proportional to the constraint itself, which vanish by the equations of motion for 
the auxilliary variable if the constraint is formulated appropriately) [lfj . From the above, we can immediately read 
off two classes of static solutions- commuting matrices, and fuzzy two spheres: 

[X h X j ]=iKc ijk X k (11) 

Where in all but name, various solutions and their energetics were discussed in the introduction ® - J7J). From the 
form of H1U|) . one might think that we've exhausted the set of possible static solutions, however a slight recasting of 
the problem will show that this is not the case. Let us perform the following change of variable: 

X + = X 1 + iX 2 , X_ = X x - iX 2 ; x\ = X_ (12) 
After which, the equations of motion become: 

X+ = l[X+,[X+,X-]] + [X 3 ,[X + ,X 3 ] + 2kX + ] (13) 

X- = ±[X-,[X-,X+]] + {X 3 ,IX_,X 3 }-2kX_} (14) 

X 3 = \[X+,[Xz,XJi\ + i[X_,[X 3 ,X+]] + k{X + ,X_] (15) 
And the Gauss law condition becomes: 

[X 3 ,X 3 ] + ±[X+,XJ\ + \{X + ,X_] = (16) 

Equations (|13|l and (|14(l are adjoints of each other, and so we are left with only two independent matrix equations 
to satisfy. The fact that the equation governing either X + or AT_ has in general a hermitian and an anti-hermitian 
part accounts for the two hermitian equations of motion that we had previously. We now try the ansatz: 

[x+,x-] = \x 3 , [x 3 ,x ± ] = ±ex ± (17) 

For which eqs. 113|) - (|15|) imply 

x+ = (2k - e - x/2)ex + (is) 

X 3 = (k~9)\X 3 (19) 
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Looking for static solutions, we immediatly spot three distinct classes of solutions- commuting matrices (A = 9 = 0), 
fuzzy spheres (A = 2k, 9 — k) and something else corresponding to A = 0, 9 = 2k. That we identify the solution 
A = 2k, 9 — k with the fuzzy sphere is readily seen: 

[X + ,X_] = 2kX 3 , [X 3 , X ± ] = ±kX ± (20) 

Which is none other than the algebra of su{2) in the Cartan-Weyl basis. It should be clear that (|12(1 constructs 
the usual raising and lowering operators if the A, are taken to be elements of su(2). The third solution we uncovered 
corresponding to A = 0, 9 = 2k, happens to describe the fuzzy cylinder^ ( see G3 a review of its construction): 

[X+,X-] = , [X 3 ,X ± ] = ±2kX ± (21) 

Very briefly, we outline its construction since it is not as commonly known as the fuzzy sphere. Consider the 
commutative cylinder, parametrized by coordinates ret and (f> e [0,27r]. Functions on the cylinder can be expanded 
in terms of the following basis: 

{e m4, Uz ; /(r,0) = ^c„(r)e^ (22) 

n 

Define x + = pe 1 ^ and x_ = pe^ 1 ^, where p is the radius of the cylinder. All smooth functions on the cylinder can 
be expressed in terms of a power series in these variables. Furthermore, we have the relation x+x- = p 2 . Being a 
symplectic manifold, the cylinder naturally admits the Poisson brackets: 

r, i dfdg dgdf 

97^-3^ (23) 

where / and g are arbitrary functions on the cylinder. The Leibniz rule implies that all Poisson brackets can be 
generated from the following elementary brackets: 

{t,x±} = ±ix± , {x + ,x-} = (24) 

The relation x + x^ — p 2 is now a central element of the Lie algebra of smooth functions on the cylinder (with the 
Poisson brackets assuming the role of the Lie bracket): 

{r, x + x-} = {x±, x+x-} — (25) 
Derivatives of functions on the cylinder can now be effected by the action of the brackets: 

dlf = {r, {r, /}} ; 3 2 f = {x + , /}} (26) 

One obtains the non-commutative cylinder by 'quantizing' the Poisson structure of the commutative cylinder. This 
is effected through the prescription {,} — > j[, ] where A is the non-commutativity parameter. Although one could 
convince themselves of the sense in this prescription through a number of heuersitic arguments, this is in fact a very 
well defined and well motivated prescription for constructing non commuting geometries |33| . With this in mind, we 
can write down the structure relations for the fuzzy cylinder as: 

[t, x±] — ±Ax± , [x + ,x_] = , x + x_ = p 2 (27) 

Thus taking A = 2k, we see that (J2U) is in fact the fuzzy cylinder if we make the identifications A3 = r and 
X± = x±. Note that the equations of motion (|21J1 do not impose any relations among the Xi other than the algebra 
(|21|l . The restriction A + A_ = p 2 is something we impose in order to make an identification with the fuzzy cylinder. 
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From the relation [X3, A±] = ±2kX±, one can immediately deduce that the spectrum of X3 is integer multiples of 
2k, and the action of X + in the basis where X 3 is diagonal is that of p times the elementary shift operator: 

X 3 \n) = 2nn\n) -> X + \n) = p\n + 1) (28) 

(X- performs the opposite shift). However we have just uncovered a very important issue concerning this solution, 
namely that any representation of the cylinder algebra will necessarily be infinite dimensional. This should be apparent 
from (|27Jl . which does not describe a semi-simple Lie algebra (re-expressing the algebra in terms of X\ and X2, we 
see that they form a proper ideal in the algebra). Hence any unitary representation (i.e. one in which the generators 
Xi are hermetian- a structural requirement of this model) will necessarily be infinite dimensional. There were two 
crucial aspects of the trace which we used in order to arrive at the equations of motion (|10|l . namely the positive 
definiteness of the trace norm, and the cyclicity of the trace. The first property was used to deduce the equations of 
motion from the statement that the action is extremized by all variations in the fields: 

Tr (lx~ SXl ) = oysx > j£ = ( 29 ) 

And the second property was used to compensate for the non-commutativity of the SXi with the Xi to bring all 
the variations to one side of the trace. It is this property that fails in general for infinite dimensional matrices (c.f. 
Tr[x,p] 7^ where x and p are a canonical pair). However one can show that if taken carefuly, the equations of motion 
Ql 01 remain unchanged in the limit of infnitc dimensional matrices, except for one important respect- one finds that 
the model now admits a 'central extension' wherein the equations of motion become |lOj: 

X t = [X,, {[X^X,] - ine l3k X k )] + cj (30) 

One finds that this central extension allows for solutions corresponding to the fuzzy plane and the warped 
fuzzy plane, but this shall not concern us here. We are interested in the case where the c$ all van- 
ish which permits the solutions that we have hitherto derived: the fuzzy sphere and the fuzzy cylinder. 

One should note that if the fuzzy cylinder and the fuzzy sphere solutions are to be meaningfully compared (or 
in the situation we're interested in, dynamically interpolated), they have to be represented in the same space (in 
our case, the usual Hilbert space 7i of square integrable functions on C |lC||l7||). However it should be clear that 
su(2) does not posses any infinite dimensional unitary irreducible representations. That is any infinite dimensional 
representation of su(2) is necessarily reducible, which means that it is going to contain an infinite number of finite 
dimensional fuzzy spheres. With this in mind, let us compare the energy densities of the fuzzy cylinder solution with 
an infinite collection of identical spin j fuzzy spheres. One finds after evaluating 10, the energy density e (= E/Trl) 
of the fuzzy cylinder is given by: 

2 

e-cyi = -k 2 P 2 T (31) 
And the energy density of the infinite collection of spin j spheres is given by: 

e sph = -^K i j(j + l)T (32) 

Furthermore, setting the radius of the cylinder to be equal to the radius of any one of the fuzzy spheres (p 2 = 
K 2 j(j + 1)) we see that the two energy densities are given by: 

2 1 

e cy i = , £s P h = -gM ; M = « j(j + l ) T o (33) 

Which strongly suggests the possibility of a transition between these two solutions, provided the energy landscape 
of this model does not obstruct the transition in any way. In and [ll| the energy landscape of finite dimensional 
solutions were carefully studied and transitions between various combinations of fuzzy spheres suggested by the 
energetics of the model were explored, with some promising results. Although the study of the energy landscape of 
this model in the case of infinite dimensional matrices is quite a formidable one, we find that a few careful considerations 
suffice to pin down an ansatz which is an ideal candidate to model the decay of the fuzzy cylinder solution into a 
collection of an infinite number of spin j fuzzy spheres. 
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III. THE ANSATZ 

As a warm up, let us consider a transition from the cylinder to a fictitious infinite dimensional irrep of the fuzzy 
sphere. We begin in this way as the treatment is relatively straight forward, and the techniques and results we uncover 
here carry over easily to a more realistic setting- the study of the decay of the cylinder into an infinite collection of 
spheres. According to the equations of motion (|f 3[l - i|15|) . the two static solutions which our matrix model admits as 
solutions in the limit of infinite dimensional matrices are: 

[X+, X_) = 2kX 3 , [X a ,X±] = ±kX ± (34) 
which describes the fuzzy 2-sphere solution, and: 

[X + ,X_] = , [X 3 ,X±] = ±2kX ± , A+A_ = p 2 (35) 

which describes the fuzzy cylinder. In the basis where A3 is taken to be diagonal, we see that for an irreducible 
representation of the sphere, A3 has eigenvalues nn where n is an integer, whereas for the cylinder, A3 has eigenvalues 
2kji (we focus on the integer representations of the sphere for simplicity). Because of the relation A + A_ = p 2 for the 
cylinder, we see that X± = px± where the x± are the elementary shift operators {x±\n) — \n± 1); A3|n) = 2nn\n) - 
see the discussion around Of course, although the action of X± on the fuzzy sphere also raises and lowers the 

state acted on, because of the relation [A + , A_] = kX 3 , this action is weighted by the Clebsch-Gordon coefficients 
(X±\n) = c±\n± l);X 3 \n) = nn\n)). If a solution is to interpolate between these two cases, it would seem that the 
most obvious form it would take would be the following: 

X+(t) = D(t)x+ , X 3 (t) = A(t)r (36) 

Where x + is the (unitary) elementary shift operator introduced above (x + x^ = x^x + = 1), and r is the operator 
representing the non-compact coordinate on the fuzzy cylinder: [r, x±] — ±2kx±, [x + ,x_] = 0. That is, we try to 
deform our cylinder solution into the spherical solution. The operators D(t) and A(t) (and all their time derivatives) 
are diagonal with real elements, as A3 only needs its spectrum to be rescaled in making the transition from the the 
cylinder to the sphere. Similarly, the action of x± also only has to be rescaled so that instead of effecting the shift 
\n) — > \n ± 1), the action of this operator becomes that of the raising or lowering operators of su(2). Furthermore, 
recalling (|13[) - l|14l) we see that it suffices to study the equation of motion for A + as the equation for A_ is the adjoint 
of the former (see the discussion following (|16fl ). Furthermore lets say that at any given time, we force our ansatz to 
satisfy the following set of relations: 

[X+(t), X_(t)] - \(t)X 3 (t) , [X 3 (t),X ± (t)} = ±0(t)X ± (t) (37) 

Which is to say, our solution will be such that it undergoes a time dependent deformation which takes it from 
the cylinder to the sphere. With this last qualification, we can straightforwardly evaluate the equations of motion 
(|13H - l|15ll . resulting in the following: 

A+ = (2k- 6- \/2)6X + (38) 



X 3 = (n-d)\X 3 (39) 

Where we suppress all time arguments from here on. One should take care in what follows in realising that 
A_(i) = X^_(t) — x-D(t), since D(t) is a diagonal operator with only real elements. The form of our ansatz Ij36|l 
trivially satisfies the Gauss law constraint 10 and (|16f) : 

[At, At] + i ([DX+, x_D] + [x_D, Dx + f) = (40) 

Where the first term in the above vanishes as both A and r are diagonal operators. The terms in the parentheses 
sum to zero, this time because both D and D are diagonal. Hence the structure of our ansatz already takes the 
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Gauss law constraint into account. We note in passing that the Gauss condition for finite N, corresponds to setting 
all of the U(N) charges to zero. One sees this from calculating the Noether charges of © under the U(N) symmetry 
(Xi — > UXiU^), which are given by: 



c a = Tr(X i [X i ,T a ]) 
= TriT^X^Xi]) 
= 

Where the T a are the generators of U(N). We used the cyclicity of the trace in arriving at the penultimate step 
above, thus we cannot make this conclusion in the case of infinite dimensional matrices. However this fact is only of 
passing interest to the problem at hand, with which we will persevere. Requiring the conditions lp?7| to be satisfied 
by our ansatz (|36|l implies the following: 



[X+, X-] = XX 3 -> D 2 - x-D 2 x+ = AAr (41) 
[X 3 ,X+}=8X+ -> DAx+2k + D(A- x-Ax+)(t -2k)x- = 9Dx+ (42) 

Where we repeatedly use the fact that the commutators of diagonal operators vanish, and the relations [r, x±] — 
±2kx± in evaluating the above. The specific form for (|42|l is arrived at by trying to get the the operators D and 
x+ to the left and right hand sides of each term respectively. We can satisfy the above conditions if the following 
equations hold true: 



D 2 - x-D 2 x+ - AAt = (43) 
2kA + (A-:e + A2;_)(t-2k) = 9 (44) 

Now realising that the spectrum of r is of the form 2ku, and the nature of the x± as shift operators, we can write 
down l|44[l for each (diagonal) matrix element: 



2n[A n + ( n - 1)(A„ - An_i)] = 9 (45) 

Further realising that X 3 — At only needs to undergo a rescaling in going from the sphere to the cylinder (X^ vl \n) = 
2nK\n), X^ ph \n) = n,K,\n)) we see that more than just being diagonal, picking A to be a time dependent multiple of 
the identity A = cul will do the job. With this in mind, the above condition simplifies greatly and implies: 

uj(t) = 6/2k ; A = ujI (46) 
Armed with this information, we can now evaluate the constraint (|43|l . which now implies: 



Dl - D 2 n+1 = n\6 (47) 

Now we know that the job of D{t) is to rescale x + such that X+(to) — px+ and X + (tf) = J + , i.e. to transform the 
cylinder X + to the spherical X + , which is the usual raising operator of angular momentum, over a given time interval. 
Let us take p to be equal to the radius of the sphere: p = Ky/ N(N + 1). Recall that the angular momentum raising 
operator J + = X\ + iX 2 (where [Xi,Xj] — ineijkXk), acts on states as J+\n) = Ky/N(N + 1) — n(n + l)|n + 1). We 
emphasize that we should suspend our disbelief at the factor N which has started to appear for the moment, as this 
is a warm up to the more realistic situation we shall consider shortly. We thus suppose a reasonable ansatz for D (t) 
would be: 



D n (t) = Ky /N(N + l)-g(t)n(n-l) (48) 

One can convince themselves of the factor n(n — 1) (as opposed to n(n + 1) which appears in the Clebsch-Gordon 
coefficients) is appropriate in the ansatz for D(t), by realising that this is the appropriate factor to left multiply 
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x + in order to scale it to the raising operator. This is easiest seen by writing out this multiplication in a matrix 
representation. However, one can also see this from the following- if: 



(n + 1| J+\n) = Ky/N(N + l)-n(n+l) , (n + l|Dar+|n) = D n+1 
Then at the time at which the transition is complete, D(tj)x+ = J+, so that 

D n+1 (t f ) = Ky /N(N + l)-n(n+l) 

so that D n (tf) = Ky / N(N + 1) — n(n — 1), hence our anstatz l)48[l. Clearly if g(to) = and gitf) = 1, then we 
satisfy the appropriate initial and final conditions: 



X +n (t ) = Ky/N(N+l) X+n , X +n {t f ) = Ky/N(N + l)-n(n + l)x +n = J +n (49) 
Hence, l|17jl implies that: 



X9 
2k 



— = 9{t) (50) 



We note that this is a highly non-trivial constraint to consistently satisfy, as in general any one parameter ansatz 
for a matrix equation has to simultaneously solve an infinite number of equations (for each matrix element), which 
we fortunately manage in our case by a judicious choice for our ansatz. We observe that (|50|l is consitent with our 
initial and final conditions: when A = 0, 6 = 2k which corresponds to the cylinder, g = 0. When A = 2k, 6 = k which 
corresponds to the sphere, g = 1 (c.f. lf3*5f and J33). Having satisfied both constraints in IpTfjl . we are ready to 
solve for the equations of motion. Before we do this, we summarize our findings: 

A(t) = e(t)/2K ; X 3 (t) = M r (51) 

2k 

D n (t) = K^N(N+l)-g(t)n(n-l) ; g[t) = (52) 
Thus our equations of motion (|38(l and (|39|l rather straighforwardly imply the following equations for our parameters: 

6 = 2{K-0)K 2 g (53) 
D n = [(2 K -9)-K 2 g]D n (54) 

We see that in general, (|54(l could be inconsistent- writing out explicitly what this equation implies: 

D n gn(n - 1)k 2 g 2 n 2 (n-l) 2 K 4 2 

D~ n = m 4D* = [(2k-0)9-k 9 ] (55) 

We see that the left hand side depends on n, whereas the right hand side does not. Consider the factor that appears 
in either term on the left: 

n(n — 1)k 2 1 . _ 

Because of this, we might think that since this model is only valid in the large N limit, both sides of Ij55(l have 
to vanish (which ensures that the n dependence of the equation is inconsequential). However this is an incorrect 
conclusion. N would truly have to be infinite if we are to even admit the fuzzy cylinder as a solution, and this calls 
into question the very consistency of what we are doing, which we remind the reader again is simply warming up for 
a more realistic treatment. We see that if this equation is to be satisfied for all finite n, then a safe way to do it would 
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be to demand that g evolve adiabatically: g 2 , g <C 1. That is, if we assume the adiabatic behaviour of our ansatz, 
the left, and consequently the right hand side of Ij55(l would have to be vanishingly small: 

(2k - 6)6 = n 2 g + e (57) 
Which we can safely take to imply the following for (|53|) : 



6 = 2(« - 6)(2k - 6)6 (58) 

Certainly we would not have to be too fussy about the adiabatic condition if we are only concerned with the 
evolution of matrix elements at finite n, if we took TV to be very large (but still finite). This is because the factors 
(15611 multiplying the time derivatives of g are highly suppressed. Expressing 9 in terms of dimensionless variables 
(6 = 6/k), the equations above become: 



6 = 2k 2 (1-6)(2-9)9 (59) 
g = 6(2-6) (60) 

From which we see that if (|59|l and (|60[1 are to be consistent with each other, and the adiabatic condition, then 
neccesarily k 2 -C 1. That is, through the relation 1)61) [I. if g is to vary slowly, then so must 9. This is only possible for 
weak coupling to the RR background, which is equivalent to the commutative limit of the model. We will have more 
to say on this further on. The equation of motion for 9 <|53[) is equivalent to motion in the potential: 



V(6) = 2k 2 (6 3 - 6 2 - 6 —) (61) 

Which we plot in Fig. 1. We note that we are only interested in the behaviour of the model for 6 ranging between 
and 2, the regions exterior to this domain having been included to illustrate the rather generic 'wine bottle' shape 
of this potential. If this model is to be taken seriously, we can immediately discern several salient features of the 
solution from this potential. Firstly, since this ansatz is only consistent in the limit k 2 -C 1, the energies of the 
spherical and cylindrical solutions l|33|) tend towards degeneracy with that of the commuting solution (a gas of DO 
branes), and hence should freely transmute into one another. From the potential, we see that we have uncovered 
the marginal instability in this limit of the commuting and cylindrical solutions (which correspond to 6 = 2 and 
respectively)- the instability is parametrized by the curvature at the maxima, which depends on k 2 , and is hence 
very slight. We note in passing that this treatment says nothing of the stability of the cylindrical solutions for finite 
values of k 2 . We also see that the spherical solution (which corresponds to 6 — 1) is still a marginally bound state. 
We can interpret the dynamics resulting from the anstaz in two ways- either as an instanton transition between a gas 
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FIG. 2: Schematic representation of the transition at intermediate times 

of DO branes and the fuzzy cylinder, (which would correspond to motion between 9 = and 9 = 2), or one could only 
follow the ansatz only up until the fuzzy sphere has been formed. Clearly unless there is some dissipative mechanism 
involved, the resulting post transition fuzzy sphere geometry cannot be static. However, if at the point 9 = 1, we 
were to drop all the restrictions (|36f) l|37|) and follow the solutions from then on, it is not inconcievable that a spherical 
solution might persist, except now as some undulating deformation of the sphere. Whatever the final outcome, were 
this admitedly unrealistic example to be taken seriously, it seems that we have modelled topological transitions using 
matrix dynamics in the large N limit. Of course any topology change would be occuring at infinity (spatially, as well 
as in terms of matrix elements) in this model, where it distinctly inconsistent to draw any firm conclusions. However 
we find that having acquainted ourselves with the nature of the problem, we only need to make a few modifications to 
this set up to model a much more realistic situation. We find that in the case where we study the decay of the cylinder 
into an infinite collection of spheres, we come across many of the features we just uncovered in this toy example but 
in more grounded circumstances, which is where we turn our attention towards now. 

IV. THE DECAY OF THE FUZZY CYLINDER: A CLASSICAL TRANSITION? 

We note that although the conditions (|37|l facilitated the subsequent calculations, they may overly restirct our 
problem. They were proposed so that at each point along the trajectory of the solution, we had a firm interpretation 
of the intermediate state as a deformation of the cylinder. Since we are proposing a solution which transmutes bewteen 
the cylinder and an infinite collection of spheres (see fig. 2), it seems reasonable to suppose that one may not have 
such an interpretation available in this case. It turns out that if we were to drop this condition at the outset, the 
equations of motion would imply them for matrix elements which are not involved (in a sense we will make precise 
later) with the topology change, whereas this will not be the case for the matrix elements which are directly involved 
in this transition. We find that dropping the conditions 137(1 does not reduce the amount of work we have to do 
in any sense, as the difficulty in constructing a solution is shifted from trying to make our ansatz consistent with 
this condition, to evaluating the equations of motion, which these conditions simplified. Before we can proceed, we 
take note of an important aspect of the fuzzy sphere solutions that we've glossed over till now. Consider the static 
equations of motion: 



[Xj, [Xi,Xj] - ine i:jk X k ] = 



(62) 
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We note that the fuzzy sphere solution given by [JQ, Xj] 



Xk admits the generalization 



X, = nJi + Y t (63) 
Where the Yi are any operators that commute with the generators of su{2) given by the Ji and themselves: 



[Y u J k ] = 0, [Y i ,Y j ] = 

In particular, if we consider the reducible solution given by Xi = 
representations, then an obvious operator that satisfies (|64ll would be given by: 



(r) 

kJ; , where the sum is over irreducible 



(64) 



Yi = 



(65) 



Where the cf are arbitrary coefficients weighting the P r \ which are the identity operators on the subspaces defined 
by the reducible representation labelled by r. These so called marginal deformations do not change the energy of the 
solutions, and correspond to center of mass dispacements of the individual fuzzy spheres that make up the solution. 
One can then easily construct a static situation where these center of mass displacements are such that we have an 
infinite tower of spin j fuzzy spheres along the X% direction. This is the solution that we will attempt to connect to 
the cylinder. To facilitate the process of trying to construct our ansatz, let us consider what the two solutions we are 
trying to dynamically relate look like in a matrix representation. When our solution describes the cylinder (at the 
initial time), the operator X+ takes on the form: 
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(66) 



Where we have partitioned the matrix into N by N(= 2j + 1) blocks. The matrix representation of X + , when it 
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describes an infinite tower (along the X3 axis) of spin j spheres is given schematically as: 



X + = K 
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Where c™ are the Clebsch- Gordon coefficients given by c" = + 1) — n(n + 1). More succinctly, we can say 

that the when X + is describing an infinite tower of spin j spheres, it is a block diagonal series of the spin j raising 
operator J + . We take the only nonzero deformation Yi to be Y3, which will have the general form l|65[l . Namely, the 
final collection of fuzzy spheres will all be displaced along the z direction. Since X3 will be a diagonal operator in the 
initial and final state, we can assume that the full time dependence of the operator will be given as before: 



X 3 (t) = A(t)T (68) 

Where r is the non-compact coordinate of the cylinder, with A also taken to be diagonal. If one were to consider 
the j x j matrix which transforms the diagonal blocks in (|66|) to those in (|67|) : 



/0 1 
1 
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0/ 



fo 4" 1 
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(69) 



We see that an ideal candidate would be the matrix d(t), similar to the one defined previously but now, a finite 
dimensional matrix (recall that if we set the radius of the cylinder to that of any of the individual spin j spheres we 
are trying to evolve to, then p = k^J j(j + 1)): 



d m (t) := Ky/j(j + l)-g(t)m(m-l) ; -j<m<j (70) 

Notice that d-j = K\/ j(j + 1)(1 — g{t)) seems to be a redundant element of our definition, as when we multiply 
the initial block matrix in (|69|l by d(t), it's —j element does not multiply any non-zero elements and hence drops 
out. However it serves the very important function of scaling the off-diagonal block elements in (|66|l to zero at the 
end of the evolution (when g = 1)- the non zero elements being in the upper off-diagonal blocks. That is, if we were 
to define the operator D{t) as the direct sum of an infinite number of these (2j + 1) x (2j + 1) block matrices, we see 
that the time dependent ansatz: 



X+(t) = D{t)x + ; D(t) = diag(...d(t),d(t),d(t)...) 



(71) 
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can scale (|66|) to (|67|) . where x+ as before is the elementary shift operator (see discussion surrounding i|36|) ). The 
equations of motion (|13fl and (|15l) now become: 

At = [(« - S1)L> 2 - a;_(K - 0)L> 2 :r + ] (72) 



D = L>[i(a; + D 2 x_ + x_D 2 x+ - 2D 2 ) + (2k - fi)Q] (73) 

fi:= [2kA + (A-:e + A3;_)(t-2k)] (74) 

Obtaining these is a straightforward excercise, where we only use the cylinder relations [r, x±] — ±2nx±, and the 
fact that the operators A and D are diagonal (recall from the discussion around l|40(l that the Gauss law constraint 
HIGH is trivially satisfied in this case) . Evaulating the above for specific forms for these scaling operators in a consistent 
way is then the remaining challenge. Before we do this, we outline a few preliminaries that will help us further on. 
Since we are aiming to transmute the cylinder solution into an infinite collection of spin j fuzzy spheres, each displaced 
along the x% axis, it should be clear that the final form of the solution should be something like: 

X 3 (t f ) = (Br{Ul + clI r ) (75) 

Where the V are the identity operators on the 2j + 1 dimensional subspaces defined by our decomposition, and the 
C3 are some as yet undetermined displacements. Consider now the following counting scheme for basis vectors in our 
Hilbert space- let any value for the index n be given by 

n = (2j + l)w + to ; w e Z , — j < to < j (76) 

That is, w labels the block we are studying and to represents the index within that block. Consider now the 
operator SI {71J, if we define A as (c.f. our ansatz in the previous section): 

. , , 6(t) u(t)w 
Zk Znn 

where 9 and u> are as yet unspecified functions of time, and w implicitly depends on n from l|76|) . It is straightforward 
to see effect of this definition on X^(t): 

Q 

X 3n (t) = A n {t)T n = —T n + LUW (78) 
ZK 

That is the effect of A is to rescale the cylinder operator, in addition to imparting a block dependent translation 
(this translation is independent of the index m within the block- see GUJO- This is precisely what we expect of an 
ansatz that will evolve to (|75|l . Thus with this definition, we can evaluate VL as: 

&>n = £l(2j+t)w+m — Q ; m -3 

= 9 + u ; to = —j 

Where the fact that w changes when we move from one block to another accounts for the difference when n is such 
that to = — j (i.e. we are considering the index corresponding to the lowest weight state within that block). Recalling 
the definition of D from (|70|l and (|71J) . it is a straight forward matter to see that for n such that m^±j, our ansatz 
satisfies (I72|l if: 

6 = 2n 2 g(n-9) (79) 
u = -(2j + l)0 (80) 
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These conditions also imply that l|72|) is satisfied When n is such that m = —j. However in the case when m = j, 
(|72|l explicitly reads: 

[(2j + l)w + j}0 + uw = (k- 6)Dl -(k-0- uj)D 2 n+1 (81) 

Which implies (recalling that the index w is one higher for the state n + 1 when n is such that m = j, that is the 
state n + 1 resides in the next highest block representation) : 

uJK 2 j(j + l)(l-g) + 2j(k - 6)n 2 g = (2j + l)w0 + u)ii + j9 (82) 
The conditions i|79|) and (|80|l imply that several cancellations occur, leaving us with the condition: 



<jK 2 j(j + l)(l-<?) = (83) 

Before we can interpret this, we use (|8U|) to solve for u> consistent with the initial conditions 8 (to) — 2k, U)(to) = 
(which come from requiring A (to) = I (I77|l 'l: 



u = -(2j + l)0 + 2/«(2j + 1) 
= (2j + l)(2«-d) 

From which we conclude that the above condition becomes: 

K 2 (2J+l)j(j + l)(2 K -9)(l~g) -0 (84) 

Which is not going to be satisfied at intermediate times in our evolution. However we will see further on that a 
limit which we will be compelled to take will ensure that this quantity is negligible. We are still left with one more 
matrix equations of motion l)73|l to satisfy. Evaluating this for n such that m ^ ±j, we get: 

D n = £>„[(2k-0)0- n 2 g] (85) 

We recognize this immediately as (15411 from our previous section. We can argue just as we did then, namely that 
since this equation implies the following: 

Dn qm(m — l)n 2 q 2 m 2 (m— 1) 2 k 4 q , , , , „ 

^ = ~ 2D 2 ~- \ D 4 = [(2n-9)0-K 2 g] ; n = (2j + l)w + m , m ± ±j (86) 

That in order for this equation to be saisfied for all m ^ ±j, our evolution should be adiabatic (g < 1, j 2 < 1), 
and the right hand side should be vanishingly small |34|. That is, just as in the previous section: 

n 2 g - (2k -6)8 = e (87) 
which we can safely take to imply for the equation of motion H79|l the following (defining 9 — 9 /2k): 



9 = 2k 2 (1-9)(2-9)9 
9 = 0(2-0) (89) 

Which are exactly the same equations as l|59^) (|5T7|) , which implied that our entire ansatz could be boiled down to the 
behaviour of a single parameter which appeared as if it was moving in the potential l|61() indicated in fig. 1. Hence it 
appears that in order for our ansatz to be self-consistent, the adiabatic condition translates into the fact that our RR 
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field strength k must be very small. However we run into an inconsistency with this demand. Examining the factors 
weighting the time derivatives of g in l|86|): 

m(m- 1)k 2 _ 1 

D 2 r i 

m Lm(m-l) ~ 9\ 

We see that we cannot safely assume that these vanish uniformly. Recall that since in the previous section, N took 
the place of j which in the (unrealistic) infinite TV limit (provided we looked at matrix elements indexed by a finite 
value of n), this factor would vanish. Since we are attempting to study the transition to finite dimensional spheres, we 
see that in general this factor will not be very small for a great many matrix elements. This would not be a problem 
if the time derivatives of g were smaller than any other scale of interest in the problem, since it is the product of the 
time derivates with these factors that we'd like to be vanishingly small. However in this problem, since the right hand 
side of i|86|l is taken to vanish almost completely (and independently) through the identification (|89p. the adiabatic 
condition implies that g, g 2 ~ k 2 (|88|)H89|1 . Hence to ignore this would be inconsistent since we're clearly considering 
terms of this order in (|88|) . However this situation becomes even worse when we consider the behaviour of the factor 
(|90|l when m — > ±j: 



m(m — 1)k 2 1 



D 2 m \^X-g] (1-.9) 

L m{m—l) ^ \ 



(91) 



which clearly diverges as g — > 1 35] . This corresponds to points in our trajectory becoming inconsistent with 
the equations of motion as we near the spherical solution. That is even if we could somehow enforce the adiabatic 
condition consistently, the approximation breaks down violently at matrix elements that are near the poles of the 
fuzzy spheres we are trying to evolve to. That is, it appears as if there is no way to ensure a consistent solution to 
(|73ll with this ansatz. Recall that what saved us from running into this problem in the case of a (fictitious) infinite 
dimensional representation of the sphere was the fact that the factor (|56[) . for finite n was extremely small in the 
large N limit. However such considerations completely ignore the matrix elements which are somehow involved in 
the topology changing process. We see that in the more realistic example of studying the decay of the fuzzy cylinder 
into an infinite collection of fuzzy spheres, our ansatz simply does not work as a solution to the classical equations 
of motion up to order k 2 . However upon closer examination, we find that this attempt yields to several interesting 
resolutions, which we address now. 

V. TOPOLOGY CHANGE: CONTACT INTERACTIONS AND QUANTUM TRANSITIONS 

In order to have a clearer train of thought in this section, we briefly summarize the findings of the last section. We 
started with the observation that the ansatz given by: 

X+(t) = D(t)x+ ; D(t) = diag(...d(t),d(t),d(t)...) 



dm(t) := Kyjj(j + 1) - g(t)m(m - 1) 

„ / s \ / \ a , s 6(t) u)(t)w 

X 3 (t) = A(t)r ; A„(0 = v i + ^"^ 

2k Zkh 

n := (2j + l)w + m 

satisfied the equations of motion (|72ll - 174|l for almost all matrix elements if the we made the identifications 

= 2K 2 g(n-6) 



w = -(2j + l)9 
g = 0(2-6) 

i > g , g 2 

This corresponded to the following equation of motion for 8 — 8 /2k: 
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9 = 2k 2 (1-9){2-9)9 

Which implied that in order for our ansatz to be consistent, k 2 <C 1- That is, our entire ansatz could be boiled 
down to the motion of a single parameter 9 in the potential (|61|) given in figure 1. However we found that this ansatz 
violated the equation of motion for X3 (I72|l for matrix elements indexed by n such that m = —j, up to a term of order 
k 3 Ij84(l . That is, this ansatz fails to solve the equations of motion for X3 at matrix elements which are the lowest 
states of a given block representation, up to order n 3 . This should not concern us too much in light of the condition 
k 2 <C 1. However we found that according to (|86[1 . in attempting to satisfy the equation of motion for X+ (|73|) . we 
could only ensure consistency with our ansatz for matrix elements such that the factor (|90|l given by: 



m(m —- 1)k 2 1 

d 2 = r ZS+S of 

happen to be vanishingly small. This arises from the realisation that although the right hand side of (|86|) does not 
depend on n, the left hand side does, unless both sides of the equation vanish (or are vanishingly small). We attempted 
this for the right hand side through the identification K 2 g — 9(2 — 9) = e, where we take e to be small independent of 
k. As for the left hand side, we require all time derivatives of g to be small, but the best we could manage consistent 
with the equations of motion for 6> <|88[1 189(1 was g 2 , g ~ k 2 , which also has to be taken to be small. However since 
we are clearly not neglecting terms of order k 2 , we need the factors multiplying g 2 and g in (|86[) to vanish. However 
these factors are not unifromly small, and for matrix elements close to the poles of the block representations we are 
considering (m — * ±j), this factor behaves as: 



m(m-l)n 2 _ 1 1_ 

Dl [JU±lL- g ] 1-g 

Which diverges as we head towards the spherical solution (g — > 1)- even within the adiabatic approximation 
we cannot satisfy the equations of motion for matrix elements that approach the highest or lowest weight states 
within a block. Hence in the limit that k 2 <C 1 our ansatz satisfies the matrix equations of motion only for 
those matrix elements of (|73|l which do not represent regions close to the poles of the spin j spheres we are 
trying to evolve to. However the matrix equation l|72|l is satisfied consistently by our anstaz for all matrix 
elements, the only violations occuring exactly at the poles and are of negligable magnitude in the above limit. 

Although our ansatz fails as a solution to the classical equations of motion, we find that the manner in which it 
fails to be rather interesting. Reconsider (JSSJ: 

Dn = Jgmjm- lW _ ^Vl)V = _ _ 2 

D n 2DI 4L>4 Y\ ) y\ V ) 

Our ansatz failed to solve this equation of motion only for the matrix elements close to the poles (m = ±j), and 
for times such that g — * 1. Since we must have g 2 , g ~ k 2 , we see that in this limit, the left hand side of the above 
diverges as l(9*Tjl: 



Suppose now we asked the question, what addition to our equations of motion would make our ansatz consistent? 
We see that in the limit k 2 — * 0, with g such that g(to) — 0, g(tf) — 1, (|94|l can be made to converge to the functional 
(if the limit k 2 — > is taken in a particular way): 
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Where A is some (finite) constant. So if instead of the equation of motion (|13|) (which gave rise to equation (|73|) ~) 
we had the equation of motion: 

X+ = ±[X+,[X+,X-]] + [X Zl [X + ,X 3 \ + 2nX + ] +XS(t-t f )MX + (96) 

where M is a diagonal matrix such that M n = 1 if n is such that m = ±j and M n = otherwise, we would have 
in addition precisely this delta function term on the right hand side of (|93|) . Our ansatz would now be a consistent 
solution to this new equation of motion. One could interpret this as arising from a contact (Yukawa type) interaction 
added to the action (0): 

S int = Jdt TriXiYijXj) (97) 

With y 3i = Y12 = Y 2 \ = 0, Yii = Y22 = XS(t—tf)M giving us our desired modification (since the equation of motion 
for X3 (|72f) is satisfied by our ansatz, we would like Y^ to vanish so as not to alter this equation). Although this might 
seem to be an ad-hoc imposition, it is not inconceivable that such an interaction might arise from more fundamental 
string theoretical considerations. Certainly considering the effect of additional background fields assuming a non-zero 
value in J^J) might effect this type of term. For instance, were the dilaton background to assume a non-zero expectation 
value, a term like the one above would be induced, except in this case, one would have: 

Yn = m\hil (98) 

Although the form of the operator M might seem a bit contrived, it is in fact far from contrived. All this extra inter- 
action term serves to do is to provide a source like term for the matrix elements which lie at the poles of the spheres, and 
as such could be considered as particles which mediate the topology change, or perturbations which cause the topology 
change. If one were to consider the action of the Yukawa interaction as a series of perturbations, we see that these per- 
turbations are what cause the 'rips' to occur. There is also a suggestive analogy between this and the action of vertex 
operators on the string worldsheet, but this analogy at present is nothing more than qualitative- our matrix model has 
very little to do with conformal field theory in this guise. Hence, the main conclusion we care to draw from this is that 
topology change (as far as this ansatz is concerned) requires our model to be perturbed somehow, either in an ad-hoc 
manner or through the mediation of some other degrees of freedom. Nevertheless, we feel that as far as uncovering the 
existence of topological dynamics in this model as a classical phenomenon, our ansatz seems to have risen to the task. 

We take note of some points that we've been somewhat cavalier about in the treatment above. Careful consideration 
of (|93|) shows that in general, if j is very large, then the factor behaves as l|92|l for many matrix elements near the 
poles of our block representation. Of course one could always modify the operator M to account for this, but this 
would make the problem start to look rather contrived. One thus concludes that the spin j spheres we are trying 
to evolve to cannot be too big. Thus we have a play-off between the limit k 2 — > and the fact that j cannot be 
too large. Whereas the former is trying to make the resultant spheres commutative, the latter is trying to keep 
the discrete structure manifest. The other fact that we've been somewhat cavalier about is that in order for the 
considerations going into (|9"5|l to be consistent (and for the violations of il9"3|) to truly occur only close to the poles), 
we really need k 2 to be very small indeed. Since the non-comutativity parameter in the structure relations is 
always k, we will still be dealing with fuzzy geometries, but only in a limit where the parameter tends towards the 
commutative limit. Note that there is no apparent paradox in our findings- topology change as we've uncovered it 
is a distictly non-commutative phenomenon. The main parameter that governs the physical non-commutativity of 
any fuzzy geometry is the dimension of the matrices being dealt with, as the algebra of N x N matrices will tend to 
the algebra of C°° functions only in the infinite N limit. Hence we can interpret the limit taken in our ansatz as 
that corresponding to a marginally commutative cylinder decaying into a series of distinctly non-commutative spheres. 

In addition to the interpretation of our ansatz as describing classical topology change in the situation where our 
model is acted on by a contact interaction with some external field, one has an interpretation of this solution in terms 
of a quantum transition readily available. Recall that our ansatz satisfied all equations of motion except those for 
(|73ll and even then, in the limit k 2 — ► 0, only for those matrix elements that corresponded to the poles of the spin j 
fuzzy spheres. So if we were to consider the statement that the action is extremized with respect to variations in the 
Xi, we see that: 



(99) 
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only for matrix elements that do not correspond to m — > ±j and for times such that t ^ tf. Let us consider the 
action evaluated for trajectories close to our ansatz: 



S = S ansatz + J dt Tr—SX, (100) 

The action for our ansatz will scale with the number of independent degrees of freedom, which for our ansatz goes 
as 3j-/V, where j is the dimension of the block partitions and N (which will tend to infinity) corresponds to the number 
of these blocks. However by 1|99|1 and the discussion above, although the first order variation vanishes for a great deal 
of matrix elements, the elements for which it does not vanish go as 2N. We arrive at this from the realisation that we 
only violated the equations of motion at the boundaries of the block representations, for times close to tf, and hence 
scales as twice the number of blocks (the factor 2 comes from accounting for the analogous equation for So if 

we were to consider the transition amplitude given by: 



JvXe--i: dtS W (101) 

where the limits on the integration correspond to path integrating over all trajectories that begin on the cylinder 
solution and end on the infinite collection of spin j spheres. One could have one of two situations to contend with- 
either a true classical solution exists which our ansatz fails to capture, or such a solution does not exist. If the prior 
case is true, then we have nothing to prove. However if such a solution does not exist, then consider the contribution 
to the path integral coming from the trajectories close to that described by our ansatz- one can conclude from the 
considerations above that: 



~0(3jN) > — \ ansatz ~0(2iV) (102) 

OA 

One recognizes in these inequalities, the conditions for the path integral to be dominated by this approximately 
stationary point. Hence in the absence of a truly classical path interpolating between the fuzzy cylinder and a 
collection of fuzzy spheres, our ansatz will be among the dominant contributions to this path integral. That is, one 
can also interpret our ansatz as the leading contribution to the quantum transition amplitude between the two static 
solutions. 



VI. CLOSING REMARKS 



Although topology changing dynamics has been a long suspected component of matrix mechanics, an explicit 
demonstration has till now proven to be somewhat elusive. Inspite of its apprerance, the existence of such a demon- 
stration is not a trivial problem to be left forgranted, as finding such a solution will answer the question of whether 
or not such topology change is a classical or quantum phenomenon. As we have hopefuly demonstrated, using matrix 
models for their ability to model the dynamics of non-commuting spaces has shown topology change to be primarily 
a classical phenomenon in the large N limit. One should interpret classical in this sense with care, as we are already 
implicitly working in a highly 'quantum mechanical' framework in proposing spacetime non-commutativity. The more 
apropriate statement being that there is no need to 'third quantize' our model to see topology change. There are 
many avenues for future investigations in light of these results. We feel that matrix models might be a highly tractable 
bootstrap into a viable theory of spacetime at the Planck length, and perhaps will imply what one already suspects 
is the case from the collective experience of string theory- that geometry is a concept derived from the physics of our 
degrees of freedom. 
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As for canonical quantum gravity, the metric tensor and its momentum conjugate with respect to the Einstein- Hilbert 
action are taken to be a (constrained) canonical pair. Although strictly speaking spacetime is still taken to be commuting, 
the metric tensor and its momentum conjugate satisfy Heisenberg inequalities, which effectively determines an energy 
dependent minimal length scale. The formation of singularities above a certain energy scale imply that there is an absolute 
minimum length scale in this theory. Slightly different considerations in semi-classical general relativity, can be shown to 
derive Heisenberg inequalities for spatial coordinates^, which imply a non-commuting structure to spacetime and hence 
imply non-commutative geometry. 
[31] See |20| for a demonstration of how T-duality as an a priori assumption in and of itself implies the appearance of a 

regularizing minimal length in field theoretical amplitudes 
[32] See and references therein- [l5f in particular. For a less mathematically cumbersome treatment, which arrives at the 
conclusion that U(N) as N — > oo should contain as subgroups, the diffeomophism groups of many 2-dimensional surfaces, 
see llal 

[33] See \\M and references therein for a lucid introduction to the quantization of classical manifolds. 

[34] Note that we need not consider this equation for m = 0, 1 as the time derivatives of D n vanish for these elements. 
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[35] We wish to point out that were we to consider the equations of motion arising from (1731 for the case when m = j, we 
arrive at the same conditions as above l|88|HI89|l . In the case when m = —j, we find that the above conditions will only 
solve the equations of motion up to a term of order k 2 - which we cannot take to vanish. Hence at m = — j in addition to 
the divergence of the factor l|91|l . we have an additional contribution to the faliure of this ansatz. 



